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Module temperature has a role in determining a PV module's performance. 
The purpose of this paper is to estimate the Joule heating in a photovoltaic 
(PV) module by comparing during PV-On (electricity generation) and PV-Off 
(without electricity generation). Joule heating was less evaluated due to 
simplifying formulation, which is easier to implement in experimental 
observation as proposed in this work. The experiment collected the 
temperature distributions of the PV module during PV-On and PV-Off. PV 
module temperature distribution follows the normal distribution curve as the 
irradiation uniformity pattern of the solar simulator has a slight <0.3 °C 
difference between PV-On and PV-Off. Joule heating slightly increased the 
PV module temperature by 0.53 K/A, proportional to the irradiances. Joule 
heating has increased almost seven times from 2.65 W at 700 W/m? to 18.07 
W at 1000 W/m’. Joule heating might slightly increase the overall thermal 
conductivity and slightly decrease the thermal resistances. It might affect the 
heat transfer. This research may improve the procedures prediction of PV or 


photovoltaic-thermal (PVT) collector temperature by considering Joule 
heating. 
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1. INTRODUCTION 

The sun is an abundant source of energy available on earth. Based on the latest research, conversion 
of solar energy into electricity can reach 24% efficiency using mono-crystalline silicon (mc-Si) of photovoltaic 
(PV) modules [1] and into thermal energy by 70% efficiency using water-based flat-plate thermal 
collectors [2]. The rest of the unused solar energy is reflected or becomes heat lost by convection. Theoretically, 
PV module performance decreases with increasing module temperature or vice versa [3], [4]. For performance 
stability, converting heat into valuable energy can control the PV module temperature. That is one of the 
reasons for combining a PV module with a thermal collector into a hybrid photovoltaic-thermal (PVT) 
collector [5]. The conversion of solar energy is up to 80-90% [6]. The performance of the PV module and 
thermal collector influence the total performance of the PVT collector [7]. 

Here, PV temperature has a role in determining optimal performance. Many researchers evaluated the 
PV module or PVT collector performance for a long time. The actual measured PV temperature is the effect of 
irradiance and Joule heating simultaneously. Joule heating is the heat generated whenever a current passes 
through a resistive material [8]. The thermal performance of a hybrid PVT collector is also likely to be affected 
by Joule heating [9], [10]. Joule heating effect exists when the hybrid PVT collector generates electrical and 
thermal energy as PVT-mode [11]. 
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PV module temperature is more determined through prediction. It used data from the manufacturer as 
environmental data (irradiation, ambient air temperature, wind speed). Moreover, the prediction data only used 
electrical parameters provided by the PV manufacturers [12]-[15]. If the prediction of PV module temperature 
is too low, then the actual efficiency of PV modules is often lower than the PV manufacturer data by standard 
test conditions (STC: 1000 W/m?, 25 °C, AM 1.5). There have been several methods and techniques in 
predicting PV temperature. Direct measurement by attaching a thermocouple to the solar cell from the 
back-sheet can accurately predict the PV temperature [16]. However, the temperature sensor insertion into a 
PV cell is quite challenging. The nominal operating cell temperature (NOCT) model is more accurate in 
predicting the PV temperature for the PVT system [17]. However, these models also depend on the accuracy 
of the environmental and operational data. Artificial intelligence (AI) techniques with artificial neural network 
(ANN) [18]-[20], modelling with various models [21]—[23] are satisfactory in estimating the PV module 
temperature. However, all these techniques have not demonstrated the physics interaction of input parameters, 
such as environmental and operational data. A model to estimate PV temperature involving wind speed effects 
gives quite reliable results [24]-[26]. However, the model is very dependent on the availability of wind speed 
data, in addition to data storage memory. Estimating the radiation and module temperature have used the 
algebra equation and non-linear parameter techniques [27]. The PV system must be on maximum power point 
tracking (MPPT). However, these methods require certainty data with quite complex calculations. In addition, 
Joule heating were not involved in the discussion. 

Joule heating has difficulty observing outdoors for a long measurement. Joule heating requires the 
ideal conditions such as; high irradiation, long duration of the steady conditions and slight temperature 
differences between inlet water and ambient air. The system must be in a steady state to be more controlled for 
getting an accurate result. While most of the indoor ones only used a small PV cell, it was not easy to detect 
and measure Joule heating. Thus, several studies have examined Joule heating by simulation rather than 
experimentation. A finite element model using COMSOL Multiphysics simulated the PV module temperature 
and Joule heating [28]. The PV module temperature was corrected using simulation for a constant current 
density. However, Joule heating was neglected for more complex simulations reasoning less contribution 
comparing irradiances and memory concerns. The comprehensive energy distribution model [29], the novel 
thermal model [30] also the electrical loss prediction of modules (ELMO), and the outdoor measurement 
method [31] investigated the electrical and thermal performance of PV modules including the series resistance 
loss from Joule heating. The model determined the amount of incident solar energy lost. However, there was 
no detailed discussion regarding Joule heating. The heat transfer model was developed based on the ambient 
condition [32]. It can implicitly involve the Joule heating in the simulation. Unfortunately, Joule heating has 
not been the main focus to examine in more detail. Finally, Joule heating did not involve in the model 
development. The 1D numerical model predicted the cell efficiency and temperature during operation [33]. 
This model encompassed every heat mechanism occurring in a solar cell. 

Research on the system performance of PV or PVT modules through module temperature is still quite 
interesting especially involving Joule heating. From previous studies, it seems that there have been some 
researches that specifically discuss the behaviour of PV or PVT module temperature as an effect of Joule 
heating. Although all studies have examined the PV module temperature, electrical and thermal efficiency; 
however, Joule heating has not been explicitly discussed [32]. The experiments could observe Joule heating 
but neglected it due to simplifying formulation. For continuing the previous study [34], the present work is to 
validate the Joule heating experimentally by investigating the thermal distribution of PV cells under PV-On or 
generating electricity and PV-Off or without generating electricity. The current study purposes in estimating 
the value of Joule heating which affected the PV module temperature during PV-On. 


2. METHOD 
2.1. Experimental setup 

Since the ideal conditions are hardly difficult to achieve outdoors, the experiments have been done 
indoors. Figure 1 shows the indoor experimental configuration using a halogen solar simulator. The halogen 
solar simulator can control the significant parameters properly. The comprehensive dimensions, designs, and 
characteristics of the simulator were previously discussed [35]. Before the experiment, the irradiance level was 
adjusted by the sliding regulator (Matsunaga, SD-1310) to the simulator and measured by a pyranometer (MS- 
42, Eko) at the midpoint of the PV module mounting surface (mc-Si, GT434 type, KIS Solar Japan). Then, the 
PV module changed the pyranometer position, connected to the simulator [35], sliding variable resistor (Berco, 
Rheostats), and analogue-digital-converter (ADC). The sliding variable resistor Rmpp adjusted the MPPT [36], 
[37]. The ADC (WEIC, Fuji electric) accumulated and converted the electrical data, such as voltage [V], 
ampere [ A ], and power [W]. The PV cell temperature Tpv(1-8) from the T-type thermocouple connected to the 
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data logger (GL220, 10 channels, Graphtec) with a 30 s sampling time. The room temperature T, as the ambient 
air temperature T, was set to 25 °C and recorded by the thermograph. 
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Computer 
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Figure 1. Indoor experimental configuration 


2.2. Detailed measurement position 

Figure 2 shows the detailed measurement position of the PV module. Figure 2(a) shows the PV cells 
coordinate at the rear surface to measure the temperature of each cell. As shown, there are four rows, each row 
with eight cells. The outer dimension of the PV module is 380x350 mm, giving the total area Apy is 
0.133 m*. The PV module had 17 cells with 35x 160 mm for each cell. Each cell gave an effective area of about 
0.0054 m°. The total area of the PV cells A,; was 0.091 m° or approximately 69% of the total area of the PV 
module A,,. Figure 2(b) describes the coordinate position of the halogen light sources. As a reference, the 
lower-left corner was the centre (0.0). The previous studies have discussed the detailed dimensions, 
construction, and characteristics of the solar simulator and the PV module [27]—[29]. Table 1 shows the specific 
measurement positions. 
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Figure 2. Detailed measurement: (a) PV cells coordinate and (b) reference position with halogen bulbs 
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Table 1. Summary coordinate position of PV cells and halogen bulbs 
Row Cells Position in PV Module Halogen Position in Solar Simulator 
1 2 3 4 5 6 7 8 1 2 3 4 
A (40;290) (75;290) (105;290) (150;290) (190;290) (225:290) (260;290) (300;290) (60;295) (145;295) (235;295) (320;295) 
B  (40;215)) (75;215) (105;215) (150;215) (190;215) (225;215) (260;215) (300;215) (60;215) (145;215) (235;215) (320;215) 
C (40;130) (75;130) (105;130) (150;130) (190;130) (225;130) (260;130) (300;130) (60;135) (145;135) (235;135) (320;135) 
D  (40;55) (75:55) (105;55) (150;55) (190;55) (225;55) (260;55) (300;55) (60;55) (145;55) (235;55) (320;55) 


2.3. Measurement and data collection 

For the study purpose, an experiment on the PV module temperature had conducted during PV-On 
and PV-Off. That is to determine the difference in temperature distribution between PV-On and PV-Off due to 
the Joule heating on the PV module. The measurement of each cell temperature took the rear surface of the PV 
module. This method allowed the light to reach the PV cell surface without being blocked by thermocouples. 
As a benchmark, measurement took 32 points. The steady-state measurements were made during MPPT with 
Rmpp™ 30; 33; 37; 41 Q for the fixed irradiances 1=1000; 900; 800; 700 W/m”, respectively. Due to the limited 
data-logger ports (only ten ports), measurements were made row by row, with the same conditions and 
treatment for each row. An experiment used eight ports to measure the temperature of each cell, the remaining 
two more ports for voltage V,,, and power P,, of the PV module. The experiment started from row A with 
I = 1000 W/m”. For measurement stability and accuracy, each endpoint of the T-thermocouple was given a 
speck of wax and affixed to the measurement point, then covered with a thin transparent tape. First, turn on the 
simulator. Then, connected the circuit for PV-On. After 30-40 minutes, the system reached the steady-state. 
With 1-minute intervals, measuring 3 data were taken for stability. Then, the circuit was cut off for PV-Off 
without turning off the simulator lights and left for about 10 minutes before repeating the same measurement. 
The same treatment and conditions repeated the experiments for rows B, C, D during PV-On and PV-Off for 
I=900; 800; 700 W/m. Table 2 shows the summaries of measurement and data collections. 


Table 2. Summary of measurement and data collections 


Row Points Activities (minutes ; 2 

EXP. (see Figure 2) (see Figure 2) 0-40 40-45 o A 60-90 Irradiances (W/m) 
1 A 8 A m=l1,2,3,...8) Steady 3data Steady 3data cooling 1000 — 900 —> 800 - 700 
2 B 8 B(@=l1,2,3,...8) Steady 3data Steady 3data cooling 1000 > 900 > 800 > 700 
3 C 8 C,(n=1,2,3,...8) Steady 3data Steady 3data cooling 1000 — 900 > 800 > 700 
4 D 8 D,@=1,2,3,...8) Steady 3data Steady 3data cooling 1000 > 900 > 800 > 700 


2.4. Theoretical and practical calculation 

As assumed previously, the effect of Joule heating resulted in a rising in temperature. The investigation 
of Joule heating needs to carry on during PV-On vs PV-Off. To illustrate the different temperature distributions 
of PV-On and PV-Off, a normal probability density function (NPDF) of each PV cell temperature Tpv,n in 
Gaussian distributions formulated as [38]: 


1 _Cpvn-n)? 1 
NPDF = f(Tpvn) = one 202 ( ) 


where; 4, 7 ando are the mean, the constant number, and the standard deviation of the PV cell temperatures 


distribution, respectively, while subscript n is cells 1 to 32. The NPDF can be proceeded using the Excel 
syntax function: NORMDIST (x, mean, standard_dev, false) [39]. 
The electrical efficiency Ne as the dependent of the PV temperature T,,,is expressed in [3], [4]: 


Ne = No|1 — 0.0045(Tpy — 25)| (2) 


where; nois the efficiency of the PV module (STC). For every 1 °C increase of T, at STC as in (2), then 77, 


will decrease by 0.0045. It is due to increased recombination losses of PV cells [29]. Practically, the electrical 
efficiency neis also given [40]; 


— Ypvlpv 


Naa (3) 
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where; Vy, Ip, and Ap,indicate as the voltage [V], the current [A] with a specific load applied [Q] to a PV 
module area [m7], respectively. The heat sources for the PV module have been identified as solar irradiance 
and Joule’s heating effect [32]. Joule heating as internal heating Qj, [W] linearly dependent on the internal 
series resistance Rsi [Q], but dependence quadratically on the current-induced I,,, [29]; 


Qin = Ri (4) 


Figure 3 shows the Joule heating conduction model from layer to layer. For simplification, assume; a). 
The system was in a steady-state with heat transfer only affected by Joule heating [41]. b). Silicon layer as a 
source of Joule heating [32], c). No Joule heating absorbed by the frame due to silicon packing, d). No loss by 
convection for Joule heating, e). PV layer temperature was the same for each point. As shown in Figure 3(a), the 
Qj, can be estimated during PV-On vs PV-Off; 


Qin = (mgiCp,gi + MevaCp,eva + Msi Cp si + MevaCp,eva + MeptCp,tpt AT py (5) 


where; m and cp are the mass [g] and specific heat [J/g.°C] of the individual layer. The subscript gl, eva, si, tpt 
are the glass, eva (ethylene vinyl acetate), silicon, and tedlar/polyester/tedlar, respectively, as shown in 
Table 3. Then, ATpy [°C] is the temperature difference between PV-On and PV-Off. 


From (5) and Figure 3, Joule heating Q;, can be formulated as (6), 


Asi Ape me Se Asi Ape aa we 


Qin = = (ksi + RevaAeva 7 an kg Agi Ta? + (ksi -> + kevafeva 7 at KepeAgi 7 Ta (6) 


where; k , A and L are the thermal conductivity [W/m.°C], the area [m°] and the thickness [m] of the 
individual layer, then the subscript gl, eva, si, tpt are the glass, eva (ethylene vinyl acetate), silicon, and 


tedlar/polyester/tedlar, respectively. To estimate the overall thermal conductivity k oy simplifying as (7); 


AT, AT, 
Apv L = or, Qin = — (7) 
pv Rov 


Qin = = 


where; Apy [m?] and Lpy [m] are the total areas and total thickness of the PV module. Then, Lyy/KpyApy is the 
overall thermal resistance Rpy [°C/W] as shown in Figure 3 (b). Table 3 shows the material properties of the 
PV module required by the formula above; thermal conductivity k [W/m.°C], specific heat c, [J/kg.°C], 
specific density p [kg/m], thickness L [m], as well as the estimated mass m [kg] and volume V [m°]. 


Silicon @-- si 
k; 
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pv (back) py (surface) | 


Figure 3. Joule heating conduction model: (a) layer to layer conduction and (b) thermal description 


Table 3. Material properties of the PV layers (summarized and estimated from [32], [40], [42]) 


Layer/Material Conductivity Specific heat Specific density Thickness Mass Volume 
(subscript) k [W/m.°C] Cp [J/kg.°C] p [kg/m] L [m] m [kg] V [mî] 

Top cover/Glass (gl) 1.04 8.40x107 2.5x10° 2.50x10% 0.404x10° 3.20x10° 
Encapsulant/EVA (eva) 0.29 2.90x10 9.6x10° 0.50x10°  0.031x10°  0.65x10? 
PV cell/Silicon (si) 142.00 6.77x10? 2.3x10° 0.24x10% 0.056x10% 0.61x10? 
Encapsulant/EVA (eva) 0.29 2.90x 108 9.6x103 0.50x107 0.031x10% — 0.65x10? 
Backsheet/TPT (tpt) 0.14 1.25x10° 1.2x10° 2.50x10° 0.404x107 3.20x10? 
Total 6.24x10% 1.906x10° 8.31x10? 
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3. RESULTS AND DISCUSSION 

Figure 4 demonstrates the PV surface temperature Tpẹ and electrical efficiency 7, from starting to 
reach a steady-state. As shown, the room temperature T, was relatively stable with an average of 25.3 °C. Each 
irradiance I had a different Typy increased, where an increase in I gave a faster Tpy rise. The fastest T,,, rise 
occurred in the collection-state (0-15 minutes) when the mass begins to collect heat, then enters the transition- 
state (15-30 minutes) starting to slow down and eventually flattens out into a steady-state (after 30 minutes). 
The effect of increasing Tp, causing a decrease in Ne (2). In the collection-state, ne decreased faster with higher 
Tpy, due to higher J. After reaching a steady-state, ne has tended to be stable. 

Figure 5 shows the bell-shaped curves normal distribution of Tp, during steady-state for PV-On and 
PV-Off. The curves are constructed based on each cell temperature with fixed irradiance J in (1), 
Figures 5(a) to 5(d) represent J=1000, 900, 800, 700 W/m?, respectively. As shown in Figures 5(a) to 5(d), the 
mean Tpy increased with increasing J. The graphs show that the difference in the temperature distribution of 
each cell decreased with decreasing I. As shown, the Tp, distribution got closer to the mean temperature. The 
normal distribution curve will be higher. However, the difference between Tpy for PV-On and PV-Off 
diminished as J decreased. 
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Figure 4. Measuring data to steady-state; PV surface temp. vs time, and electrical efficiency vs time 
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Figure 5. Distribution of Ty: (a) [=1000 W/m?, (b) [=900 W/m?, (c) 1=800 W/m?, and (d) [=700 W/m? 
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Figure 6 shows the thermal and electrical performance during steady-state. The relationship between 
PV surface temperature T,,, and irradiance 7 under PV-On and PV-Off conditions is given in Figure 6(a). As 
shown, the Typy increased with increasing J or vice versa. The difference between T;,, for PV-On and PV-Off 
was greater with increasing irradiance, and close to zero for 1~620 W/m’. Figure 6(b) shows the relationship 
between electrical efficiency ne for PV-On with PV module temperature T,,,,. As shown, the ne decreased with 
increasing T,,,, or vice versa. From the graph, the temperature efficiency coefficient u=-0.0049/°C is quite 
reasonable compared to the previous studies [3], [42], [43]. 

The complete results for PV surface temperature Tpy, standard deviation o , root mean square error € 
and Joule heating Q;, are summarized in Table 4. As shown, the irradiance decreased so that Tyv, O, €, Qip also 
decreased. Table 5 compares the present and previous results [32], such as thermal resistance R,, and 
temperature rise due to current (Joule heating) ATpy/Ipy. The lower thermal resistance Rpp in the current result 
compared to the previous result was probably due to the different PV module construction. However, the 
temperature rise due to the current AT,,/I,, was not much different. 
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Figure 6. PV performance: (a) PV surface temp. vs irradiance, and (b) elect. efficiency vs PV temp. 


Table 4. Summary results of PV-On vs PV-Off 
I= 1000 W/m? 1=900 W/m? I= 800 W/m? I= 700 W/m? 
PV-On PV-Off = PV-On PV-Off = PV-On PV-Off = PV-On PV-Off 
PV surface temperature Toy PC] 39.43 39.17 37.52 37.33 35.48 35.37 33.44 33.41 


Items 


Std. Deviation o [°C] 0.86 0.89 0.77 0.81 0.74 0.77 0.70 0.71 
Root Mean Sq. error E [°C] 0.84 0.88 0.76 0.79 0.73 0.76 0.69 0.70 
Joule heating (calculated) Qin [J] 842.52 - 744.78 - 515.72 - 123.57 - 


Table 5. Comparison to the closely related previous result [29], [32] 
Result Comparison 


asia Proposed Ying [32] 
Thermal resistance in (8), Rpy m’.K/W 0.0049 0.004 1-0.0043 
Temp. rise due to current (Joule heating) AT,y/Ipv K/A 0.53 0.54 


From the results with the condition of PV-On and PV-Off, such as; normal distribution, mean of Tyv, 
ATpv, thermal and electrical performance, and so on, three things can be discussed. First, related to the thermal 
distribution from Figure 5. The reason why the distribution of PV module temperature follows the normal curve 
is that follows the irradiance uniformity pattern of the simulator [35]. The detailed overlap coordinate position 
between halogen light sources and the PV module has shown in Figure 2. There was a slight <0.3 °C of the 
temperature difference ATpy between PV-On and PV-Off. The small ATpy may be caused by the small Lpy 
generated by the small Apy. For =1000 W/m’, the difference T,,, in the middle and at the edge was 2 °C. The 
temperature rise due to the current ATpy/lpy was 0.53 K/A. Those are close to the glass to back sheet (GB) 
module [32]. Implicitly, it might be due to irradiance as well as Joule heating [9], [10]. 

Second, related to the thermal and electrical performance of the PV module as shown in Figure 6. The 
temperature uniformity of PV cells is one of the main factors affecting electrical efficiency. With a more 
uniform temperature distribution, the electrical efficiency will be higher (2). The output power will be even 
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greater. Conversely, the PV module temperature increases, the power output will decrease. As a result, the fill 
factor (FF) decreased, which mean the electrical efficiency of PV also decreases [44]. 

Third, related to the overall thermal conductivity kp, and thermal resistance Rpy. Referring to 
Figures 3(a) and 3(b) at steady-state and constant ambient air temperature Tq, the k,, changes with load 
variations Rmpp. With fixed Apy and Lpy, when I increases, causing AT;, to increase, then Qj, also increases. 
In this case, the Joule heating might slightly increase kp, and also slightly decrease Rp, (from glass cover to 
back sheet) which affect the heat transfer [32]. The difference in operational temperature of each PV cell results 
in different performance. The PV module accuracy decreases from low-temperature toward the high- 
temperature range [45]. 

Finally, the effect of Joule heating has been proved and estimated, with the results positively 
correlated to the irradiances increase from 700-1000 W/m. For more accurate results, further experiments use 
more precise and sophisticated instrumentation. This study tries to estimate the effect of Joule heating in a PV 
module as a reference analysis for PV and PVT performances. The current study contributes to improving the 
procedure of predicting PV module temperature by considering the Joule heating. In practical application, Joule 
heating optimises the electrical efficiency of PV modules. Joule heating also increases the thermal and electrical 
performance of the PVT system as additional heat [9], [10], [15]. By involving the temperature distribution of 
each PV cell, the performance prediction and calculation of PV or PVT module temperature will be more 
accurate. For future research, the performance of a large PV or PVT system should consider Joule heating as a 
complement by improving the performance modelling tools, system design and development [46]-[48]. 


4. CONCLUSION 

The existence of Joule heating has been investigated indoors by comparing the PV surface temperature 
of each cell with the higher irradiances of 700-1000 W/m? during PV-On and PV-Off in steady-state. The PV 
module temperature distribution follows the curve as the irradiation uniformity pattern of the solar simulator 
with a slight <0.3 °C difference between PV-On and PV-Off. Joule heating slightly increased the PV module 
temperature of 0.53 °C/A, which is proportional to the irradiances. Joule heating has increased almost seven 
times from 2.65 W at 700 W/m? to 18.07 W at 1000 W/m’. Joule heating might slightly increase the overall 
thermal conductivity and also slightly decrease the thermal resistance. Joule heating affects heat transfer. For 
more accurate results, further experiments use more precise and sophisticated instrumentation. This research 
may improve the procedure in predicting the PV module temperature by considering the Joule heating during 
the PV-On and also PVT-mode. Joule heating increases the electrical efficiency of the PV module and the 
thermal efficiency of the PVT collector. For future research, consider Joule heating in the thermal and electrical 
performance of a large PV module or PVT collector by improving the modelling and control systems. 
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